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A property common to several different astrophysical sources of high-energy gamma-rays is the presence of bulk 
motion at relativistic speed. The intrinsic spectra of the nonthermal radiating particles also show interesting 
similarities, with a pronounced hardening towards lower energies. This suggests two distinct acceleration mech¬ 
anisms could be at work in these sources. At high energies, the stochastic first-order Fermi process at shocks 
seems to provide a reasonable explanation. I will briefly review the status of this mechanism before discussing 
the possibility that, at lower energies, non-stochastic acceleration in the induced electric field of a relativistic 
current sheet plays a role. 


1. INTRODUCTION 

With only one or two exceptions, the identification 
of relativistic bulk motion in an astrophysical source 
rests on the interpretation of its nonthermal radiation 
normally a featureless continuum, presumably pro¬ 
duced by relativistic electrons as either synchrotron 
radiation or by the inverse Compton scattering mech¬ 
anism. This makes it very difficult to constrain the 
underlying physics. In most cases, cut-off or break 
frequencies are only poorly defined by the data and 
the only quantity that can be used to constrain mod¬ 
els is a spectral index. Even this is not always easy to 
interpret: to be meaningful, it requires a power-law 
spectrum extending over at least a couple of decades 
in frequency. Nevertheless, the evidence accumulated 
from the spectra of shocked pulsar winds, gamma-ray 
bursts, and blazars indicates that the distribution of 
radiating particles is relatively soft at high energy 
(s = —din//din7>4, where / is the phase-space 
density of the injected particles) and hard at low fre¬ 
quencies (s<3.6). The first-order Fermi mechanism 
operating at a relativistic shock predicts indices not 
too different from those seen at high energies, but it 
fails to provide an explanation for the hard low energy 
spectra. A promising candidate mechanism in this 
regime appears to be acceleration at current sheets. 
In this paper I briefly discuss the observational evi¬ 
dence, review the predictions of the first-order Fermi 
mechanism and summarise recent ideas on relativistic 
current sheets. 


2. OBSERVATIONS 

Plerions — fille d-cen tre supernova remnants 
IWeiler and Panagialll978l] — are thought to contain 
a pulsar that powers the nonthermal nebular emission 
within the supernova bow-shock via a magnetised, 
relativistic wind. Although data on the spectra of 
sever al examples have been modelled fe.g.. IClrevalierl 
m, by far the best observed member of this class 
is the Crab Nebula. Optical and X-ray images of 
this object show features that apparently move with 


mildly relativistic velocities ~ 0.5c iHester et all 
120021 . However, an analysis of the dynamics of 
the pulsar wind responsible for the powering of 
the Nebula indicates that it is highly relativistic 
(bulk Lorentz factor P ~ 10 3 to 10 s ) and ma y also 
be strongly magnetise d iKenn el and Coronitil 11984 
Kundt^nd^Krotschecld 1980 .. j 
200ll iLvubarskvl 120031 iR.ees and Gunnlll974l| Over 
the entire X-ray energy range 2 x 10 16 to 6 x 10 20 Hz 
its spectrum is well-described by a po wer-law of pho¬ 
ton i ndex q = — din N 1 /dh\<lv = 2.11 [Massaro et all 
12000 1. In the gamma-ray region, this synchrotron 
spectrum cuts off at around 10 22 Hz, whereas to¬ 
wards lower frequencies it hardens substantially, 
reaching a photon inde x of q = 0.26 in the radio 
range (10 s to 10 lo Hz), iBaars and Hartsuiikerill9'72 ( 

iBietenholz and Kronberslll992l. 

In the case of gamma-ray bursts, there is a gen¬ 
eral consensus that a highly relativistic flow with large 
Lorentz factor F > 100 is present. Nevertheless, it is 
difficult to constrain the physics of the outflow with 
the available data. Several break frequencies are pre¬ 
dicted, and have been i dentified in the observations 
[for a review see IPira 7] 1200,511 . But direct observa¬ 
tion of a power-law spectrum over a substantial fre¬ 
quency range is elusive. The “measured” indices fre¬ 
quently represent interpolations between observation 
bands. Light curves are in some cases quite accurately 
known, but can only be used to constrain theories of 
the acceleration mechanism when combined with ad¬ 
ditional assumptions about the dynamics of the ex¬ 
pansion. Nevertheless, all models currently discussed 
conform to the trend of a hard electron spectrum at 
low energy (usually manifested as a low energy cut¬ 
off) and a softer one at high energy, although in some 
bursts, even the X-ray synchrotron emission may lie 
in the “l o w energy” regime according to this definition 
iMoran et aZ.ll2005j . 

A third class of object thought to harbour a rel¬ 
ativistic flow is that of blazars. Relativistic bulk 
speeds are required both to explain apparently super¬ 
luminal motion in their jets and by the rapid vari¬ 
ability in their gamma-ray spectra. As in the case 
of gamma-ray bursts, the spectra do not show clean 
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power-law behaviour over several decades of photon 
energy but demand a detailed analysis before they 
can be interpreted in terms of acceleration theory 
|Sambruna et al\\\ 99 fit iTavecchio et, oi.ll2002j . At ra¬ 
dio frequencies, the spectra are systematically harder 
than at high frequencies. This could either reflect the 
intrin sic particle distribution jMastichiadis and Kirk! 
119971 1 or be caused by internal absorption. Several 
individual multi-frequency flares have been success¬ 
fully modelled by using either a very hard injec¬ 
tion spectrum s ~ 3.3 combined wi th a high fre¬ 
quency cut-off |Konopelko et all 120031 or a low fre- 
que ncy cut- off and softer s ~ 4.2 high-energy power- 
law iKrawczvnski et alll2002i . 

The study of Blazar spectra illustrates the diffi¬ 
culties involved in trying to extract the intrinsic in¬ 
jected particle spectrum from an observed photon 
spectrum that has been shaped by various loss pro¬ 
cesses and by internal and possibly external absorp¬ 
tion. An interesting alternative approach to this prob¬ 
lem h as recen tly be en ad opted for FR I radio galax¬ 
ies bv lYoung et al\ (2005$ . who identify s = 4.1 as a 
characteristic power-law index of the acceleration pro¬ 
cess. This analysis provides no evidence that the low 
energy particle distribution hardens, but cannot be 
interpreted as evidence against such a trend. 



Figure 1: The power-law index s and the compression 
ratio of relativistic shocks as a function of the spatial 
component of the four-velocity of the upstream plasma 
into the shock Tm. Two cases are shown: (i) 
non-magnetised with the full Synge/Jiittner equation of 
state (solid line) and (ii) unmagnetised upstream (dotted 
line) but with an oscillating field component generated to 
the level of o+ = 1/100 (see text). 


3. RELATIVISTIC SHOCKS 


The first-order Fermi mechanism of particle accel¬ 
eration at relativistic shocks rests on the assump¬ 
tion that energetic, charged particles are transported 
stochastically through the bac kground plasma around 
a shock front [for a review see Kirk and Duffvl ll999l. 
Particles cannot be assumed to “diffuse” in space. 
This is because the diffusion equation can be derived 
only for particles whose angular distribution func¬ 
tion is approximately isotropic in the local plasma 
rest frame — an impossible requirement if the rela¬ 
tive speed of the upstream and downstream plasmas 
is comparable to the speed of the particles themselves 
iKirk and Schneideilll987l| . as, for example, at a rela¬ 
tivistic shock. 

Nevertheless, provided the particle transport pro¬ 
cess is stochastic and does not introduce a character¬ 
istic momentum scale into the problem, one can still 
expect the acceleration process to produce a power- 
law spectrum in particle energy, at least for ultra- 
relativistic particles, whose velocity can be considered 
to be independent of energy. A model of the trans¬ 
port process has to be adopted in order to find the 
power-law index, but the result appears to be rather 
insensitive to the particular choice. The kinematic 
problem of particle acceleration at a relativistic shock, 
i.e., that of finding the distribution of a collection of 
test particles undergoing small-angle, random, elas¬ 
tic (in the plasma frame) deflections in the vicinity 


of a discontinuity in the (relativistic) plasma velocity 
is well-understood. An analytic method based on an 
eigenfunction decomposition is available which gives 
the spectrum and angular dependence of the distribu¬ 
tion function at energies we ll ab o ve those of injection 
for arbitrary shock speeds |Kirk et aZll2000l . In ad- 
dition, Monte-Carlo simulations have been performed 

jAchterberg et all200ltlBednarz and Ostrowskill998j 

finding good agreement with the analytic results. 
These are illustrated in Fig. which shows the com¬ 
pression ratio and the high-energy power-law index s 
as a function of the spatial component of the 4-speed 
Tu of the upstream plasma, where F = (1 — it 2 ) -1 / 2 . 
An interesting aspect of these results is that the 
power-law index tends asymptotically to the value 
s « 4.23 for large shock Lorentz factors (or, equiva¬ 
lently, upstream Lorentz factors), independent of the 
equation of state of the plasma. This asymptotic value 
is essentially fixed by the compression ratio of the 
shock and depends only weakly on the form of the 
scattering operator used to describe the small-angle 
deflections. 

The eigenfunction expansion method enables the 
full angle dependence of the distribution to be ex¬ 
tracted, giving additional insight into the kinematics 
of the acceleration process. Both upstream and at 
the shock front itself the angular dependence is well- 
approximated by the simple expression 

/ oc (1 - p s u)~ s exp f- 1 + ^ (1) 
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Figure 2: Contour plot of the approximate angular 
dependence of accelerated particles upstream of the 
shock front (Eq.Qfor s = 4.23), plotted in the rest frame 
of the shock. v x and v y specify the particle velocity, and 
u is the three-speed of the background plasma flow, 
(which is in the positive x direction in this frame) in 
units of c. Since v x + v y = 1 for ultra-relativistic 
particles, the angular distribution for given u lies on a 
circle of radius u (shown as dashed lines for u = 0.2, 0.5 
and 0.8). The contours are linear, starting at / = 0.1 
(light) and spaced by 0.1 up to / = .8 (dark). 


where p s is the cosine of the angle between the shock 
normal and the particle velocity, measured in the 
frame in which the shock is at rest and the upstream 
plasma flows along the shock normal at speed cu. 

This function is illustrated in Fig. [3 for s = 4.23. 
Contours of / are shown as a function of a combi¬ 
nation of the particle velocity components v x and v y 
and the speed of the upstream plasma into the shock 
u. The shock normal lies along the x axis, along which 
the plasma flows. Particles with v x > 0 are stream¬ 
ing back towards the shock front; the distribution is 
rotationally symmetric about the x axis. The further 
one goes away from the shock into the upstream re¬ 
gion the better / approximates the full distribution. 
For ultra-relativistic shocks, / is a good approxima¬ 
tion even at the shock front itself, in which case those 
particles with v x > 0 are crossing from upstream to 
downstream and those with v x < 0 vice versa. For 
a given upstream plasma speed u, the angular distri¬ 
bution of the accelerated particles is determined by 
the contours intersecting the circle of radius u, cen¬ 
tred on v x = v y = 0. At low shock speeds, the 
distribution function is more or less isotropic, with 
a slight forwards/backwards asymmetry. However, 
above roughly u = 0.5, a pronounced cone emerges, 
directed back towards the shock front. In the ultra- 
relativistic limit, this cone has an opening (half) angle 


0 C given by 

0 C = arccos ( 2 ) 

It is independent of the details of the upstream trans¬ 
port, provided these can be described as a diffusion in 
angle. The details of the downstream transport enter 
only in that they have a slight influence on the power- 
law index s. The cone arises because of two competing 
physical effects. Firstly, those particles that cross the 
shock with velocity inside the cone: v x > cos 0 C are de¬ 
pleted because they have a high probability of escape 
downstream. Secondly, the angle at which a particle 
propagates in the upstream plasma reflects the energy 
and, to some degree, also the angle at which it last 
emerged from the downstream plasma. Figure| 2 ]shows 
the upstream distribution at fixed energy in the shock 
rest frame. In this case, the higher v x , the greater the 
energy gain has been since the last crossing. A large 
energy gain also implies that the potential source par¬ 
ticles are more numerous. Consequently, this effect 
leads to a steady increase in the distribution function 
as v x increases. 

Using a Monte-Carlo approach, it is possible to in¬ 
vestigate more general forms of the scattering opera¬ 
tor, whilst retainin g the effect of a non-vanishing aver¬ 
age m a gnetic fiel d [Achte rberg et aZjl200 li lOstrowskil 
1 19931 IVirt.anen and Vainiol 12005)1/ Provided the tur¬ 
bulence remains strong, little difference is found. 
However, as expected, the acceleration mechanism 
becomes less effective a s the turbulence diminishes 
lOstrowski and Bednar^ I2002I1 . because the regu¬ 
lar component of the field in the downstream re¬ 
gion quickly dr ags particles aw ay from the shock 
front iBegelman and Kirkl Il99fll . Explicit calcula¬ 
tions of particle motion in a completely random mag¬ 
netic field (with va nishing average co mpo nent) have 
been performed b y iBallard and Heavensl ^992] and 
ICasse et aQ 12002 1 . They have been used to compute 
the acceleration around a relativistic shoc k for Lorentz 
factors F < 5 [Ballard a nd Heavenslll992| and, mor e 
recently, for F < 100 tLeinoiiie^Lnd^elletieil 120031 . 
The latter find good agreement with the analytic re¬ 
sult on the asymptotic power-law index. 

Although particle transport in astrophysical plas¬ 
mas is usually dominated by interaction with fluctua¬ 
tions in the electromagnetic field produced collectively 
by the background plasma, there there are strong in¬ 
dications that two-body collisional processes (includ¬ 
ing those with the photon gas e.g., photo-pion pro¬ 
duction and Compton scattering) may be important 
for the acceleration and/or the thermalisation of en¬ 
ergetic parti cles in the inner parts of a GRB fireball 
r < 10 16 cm |Derishev et cl.lEfifKl ISternll2C)03[ . These 
generically produce much harder spectra, principally 
because an energetic particle occasionally takes on an 
uncharged “identity” as a photon or neutron. This 
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facilitates flights deep into the upstream region, en¬ 
abling the particle to profit from an energy boost of a 
factor r 2 , which is normally available to charged par¬ 
ticles only on their initial shock encounter. However, a 
hard spectrum enhances the modification of the shock 
front by the accelerated particles. This implies that 
fully nonlinear calculations will be required to assess 
the importance of collisional processes. 

On the other hand, the nonlinear modification of 
a collisionless relativistic shock does not affect the 
asymptotic power-law index. There are two reasons for 
this: firstly, isotropised, accelerated particles behave 
like a relativistic gas with adiabatic index 4/3, so that 
the overall compression ratio of an ultra-relativistic 
shock front remains 3, even when a significant part of 
the overall energy and momentum flux is carried by 
these particles. Secondly, the asymptotic power-law 
index in the test-particle picture is soft (i.e., s > 4). 
This means that it is possible to consider a Lorentz 
factor above which the test-particle approximation is 
valid, because the energy density in the remaining ac¬ 
celerated particles is indeed small. Nevertheless, a 
strong nonlinear effect can be exerted by particles of 
lower energies, whose mean free path to scattering is 
comparable to t he size of internal s tructures in the 
shock transition lEllison and Ponblell2002l . 

In order to understand how a collisionless, relativis¬ 
tic shock can form, it is necessary to identify a suit¬ 
able instability which can lead to dissipation in the 
nonlinear regime. Currently, the most promising ap¬ 
proach to this problem considers the nonlinear d evel ¬ 
opment of the Weibel instabi lity jMedvedev and Loebl 
119991. [Yang et aZ.liil994 . 19931 . which generates down¬ 
stream magnetic field perpendicular to the stream¬ 
ing motion of the plasma i.e., in the plane of the in¬ 
cipient shock. Particle-in-cell simulations of this sit¬ 
uation have been per f ormed l.laro schek et al\ [2005, 
iNishikawa et all 120031 I Silva et al 3 120031 suggesting 
that magnetic field can be generated with a strength 
<t+ of a few percent. (The magnetisation parame¬ 
ter <t+ is defined as the ratio of the magnetic energy 
density to twice the total enthalpy density (includ¬ 
ing rest mass) as measured in the downstream plasma 
rest frame). This is encouraging, since it is roughly 
the level i mplied b y spectral m odelling of GRB after¬ 
glows IPanaitescu and KumaiHEoF)^ . However, it has 
so far not been possible to identify par ticles that par¬ 
take in the first order Fermi process jHededal et al\ 
l2004ll , nor is the ultimate fate of the g enerated field 
fully understood |Medvedev et oiJl2005j . 

The manner in which magnetic field is generated 
at the shock has a strong influence on the spectrum 
of accelerated particles. However, if we are interested 
only in high energy particles of long mean free path, 
the complex aspects of the problem can be by-passed: 
the power-law index predicted by the first-order Fermi 
mechanism can be calculated simply by modifying the 
shock jump conditions to account for the generated 


field. To do this, consider time-averaged conditions, so 
that linear functions of the oscillating electromagnetic 
field vanish. The stress-energy tensor in the plasma 
frame is 


rj-1/diy 




9 


fllS 


BUB” 

4-7T 


( 3 ) 


(for notation see lKirk and Duffvl ^)99]) and the last 
term on the right hand side does not contribute to the 
fluxes across the shock front if the magnetic field lies 
in the shock plane. As a result, the jump conditions 
are the same as those of an unnragnetised fluid, pro¬ 
vided the magnetic enthalpy density B 2 /Air and pres¬ 
sure B 2 /8it are taken into account [Lvubarskvjl20031 . 
For a relativistic gas, this gives an effective adiabatic 
index 


7eff 


4(1 + 0 -+) 

(3 + <r+) 


( 4 ) 


leading to an asymptotic compression ratio of 
1 / (7eff — 1) and a relative speed of the up¬ 
stream medium with respect to the downstream 
medium corresponding to the Lorentz factor T re i = 
T-v/(2 - 7eff)/7eff- As cr+ increases, the compres¬ 
sion ratio of the shock decreases and the high-energy 
power-law softens, as shown in Fig. ^ If magnetic 
field amplification indeed saturates at <r+ ~ 1%, the 
asymptotic spectral index still remains close to 4.2. 


4. RELATIVISTIC CURRENT SHEETS 

The first order Fermi process at a collisionless, rela¬ 
tivistic shock does not appear to produce spectra with 
s $ 4, so that an additional mechanism is required in 
many sources. Since this mechanism tends to manifest 
itself at lower energies, it could also play the role of 
injecting particles into the Fermi I process. Possible 
candi date processes include acceleration by a velocity 
shear |R.ieger and Duffvll2004 IStawarz and Ostrowskil 
2002 j, the maser mechanism of lHoshino**eT^Zr [l992j 
and the destructio n of magnetic flux in the shock front 
i Lvubarskvl l2003t , as well as the second order Fermi 

f >r ocess o f acceleratio n by a turbulent wave spectrum 
Virtanen and Vainio 2005 1. Another pos sibility is ac - 
celeration at relativistic current sheets |Kir H 1200411 . 
As well as its obvious importance in the process of 
flux destruction, this possibility is particularly attrac¬ 
tive in view of the fact that field reversals of short 
length scale can be g enerated at relativistic shocks 
[Medvedev et af.ll2005l| . 

The current sheets at which reconnection and par¬ 
ticle acceleration takes place in astrophysics are rel¬ 
ativistic in two senses: Firstly, the magnetisation 
parameter, a is large and the Alfven speed v\ = 
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cy / cr/(l + a) is close to c. Secondly, the geometry 
of the current sheet at which magnetic energy is dis¬ 
sipated and, hence, the field configuration, is dictated 
by a highly relativistic plasma flow. Particle accel¬ 
eration depends crucially on both the magnetisation 
parameter and the field configuration. 

The relativistic effects associated with a large mag¬ 
netisation parameter are readily appreciated. On the 
other hand, the geometrical effects of a relativistic flow 
are more subtle. The situation is closely analogous 
to that of MHD shock fronts, which can be classi¬ 
fied into “subluminal” and “superluminal” according 
to whether the speed of the intersection point of the 
magne tic field an d the sh ock fr ont is less o r greater 
than c iBegelman and Kirkil99fllDrurvil983| . In each 
case, a Lorentz transformation enables the shock to be 
viewed from a reference frame in which it has a par¬ 
ticularly simple configuration: either a de Hofmann- 
Teller frame with zero electric field, or a frame in 
which the magnetic field is exactly perpendicular to 
the shock normal. In the case of a current sheet, the 
speed of the intersection point of the magnetic field 
lines and the sheet centre-line is important. If it is 
subluminal, a transformation to a de Hofmann-Teller 
frame is possible, leading to th e sta ndard configura¬ 
tion for a n onrelativistic sheet jBiichner and Zelenvii 
ll 989tlChenll1 992l . Alternatively, for superluminal mo¬ 
tion of the intersection point, which should be the rule 
for sheets in relativistic flows, a frame can be found 
in which the sheet is a true neutral sheet with no field 
lines linking through it. T his is, in fact, the original 
configuration considered bv ISoeise 3 fl96al . For a rel¬ 
ativistic sheet, however, it is the generic case, rather 
than a very special singular one. 

Most discussions of reconnection treat a Sweet- 
Parker or Petschek configuration in which the length 
of the current sheet in the average field direction deter¬ 
mines the dissipation rate. This is also true for recent 
analytic treatments that are relativistic i n the sense 
that the effects of large a are included iLyu barskj 
20051 iLvutikovI 1200,11 iLvutikov and Uzdenskvl l2002j . 
But the vanishing of B z in the generic relativistic case 
has important implications, since it is the linking field 
that can eject particles from the sheet, making it cru¬ 
cial for the determination of the spectrum of acceler¬ 
ated partic les, and, especially, the maximum perm it- 
ted energy jLarrabee et «/.ll2003l iLitvinenkoll 1 99fl . 


Relativistic current sheets, can extend over large 
distances along the field, depending on the nature 
of the boundary conditions. An example, d rawn 


from 

the case of a striped pulsar wind Coroniti 

199C 

Kirk and Skiaeraasen 2003, iLvubarskv and Kirk 

2001 

, is shown in Fig. |3J If we assume that recon- 


nection leads on average to a stationary field config¬ 
uration, then as the spiral pattern moves outwards, 
the linking field lines shown in the inset must move 
through the plasma at a speed sufficient to keep their 
average distance from the star constant. The striped 


Figure 3: The striped pattern of a pulsar wind. A 
magnetic dipole embedded in the star at an oblique angle 
to the rotation axis introduces field lines of both 
polarities into the equatorial plane. The current sheet 
separating these regions is shown. In the inset, an almost 
planar portion of this sheet (dashed line) is shown, 
together with the magnetic field lines, assuming they 
undergo reconnection. 


z 

T x 



spiral pattern depicted in the figure is expected to 
be established well outside the light cylinder, defined 
to be at radius r = Tl, where the corotation speed 
reaches c. In this case, the magnetic chevrons, which 
must move a distance 2nr in each rotation of the spi¬ 
ral pattern, have a superluminal speed equal to cr/r l- 
Transformation to the frame in which the sheet is a 
true neutral sheet involves a small boost in the x di¬ 
rection, and the resulting configuration has a typical 
dimension in the azimuthal direction of ~ 2nr. 

Particle acceleration in current sheets with finite 
linking field (BA has been extensively investigated 
ISvrovatskil 1 19811 . But in the generic, relativistic, 
configuration, the linking field can play no role in 
ejecting particles from the sheet. Instead, acceleration 
is controlled by the finite extent of the sheet in lati¬ 
tude, i.e., in the direction parallel to the electric field 
( E y ). This is limited not by the boundary conditions, 
but by local parameter values, as first described by 
lAlfvenl 119681 . Assuming the plasma consists of cold 
electrons and positrons, and that a 1, the maxi¬ 
mum Lorentz factor 7 max after acceleration is, 

7max = 2(7, (5) 

whereas a cold electron-proton plasma gives 

7max ~ cr for protons (6) 

7max ~ crM/m for electrons, (7) 

iKirkl 12004 1 with M and m the proton and electron 
masses, respectively. It is interesting to note that 
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in a plasma in which the magnetic field and par¬ 
ticle rest mass are in rough equipartition (er ~ 1), 
the upper limit given by Ea. (171) coincides with that 
quoted bv lLesch and Bir3 |1997| . However, this situ¬ 
ation arises only in relativistic plasmas. In the inter¬ 
stellar medium, for example, a ss 10~ 9 or smaller, in 
which case the upper limit on the energy gain reduces 
to Mv\. Standard estimates of the interstellar mag¬ 
netic field and particle density (1 pG, 1 proton/cm 3 ) 
imply that electrons can be accelerated, at most, to 
only mildly relativistic energies. In this case, and in 
solar system applications, direct acceleration by the 
DC field may be masked by particle acceleration in 
the turbu lence fed by reconnection or the associated 
shocks |Cargillll20m| . 

The picture sketched above applies only 
to quasi-steady cu rrent s heets . However, 

Zelenvi and Krasnosel’skikbl ;1 979j have shown 
that relativistic current sheets are unstable to the 
growth of the tearing mode and other i nstabilities 
are a lso likely to operate (see, for example. iDaughtonl 
1999 1 ). On scale lengths comparable to the sheet 
thickness an unsteady, oscillating component of B z 
may be generated. Thus, locally, the nonrelativistic 
picture may be relevant to the micro-structure of 
the sheet, although not in its standard 2-dimensional 
stationary incarnations. 

Particle-in-cell simulations can provide valuable in¬ 
sight here, provided they account for relativistic par- 
tic le motion. Such s imula tions been performed in 3D 
kJa roschek et !S0 129941 and IZenit,a,ni and Hoshind 
l290a| . who noted the growth of corrugations in the 
current sheet in the direction of the electric field E y 
and identifi ed them as due to the r elativistic drift kink 
instability. I.Taroschek et al\ |2004j found a very hard 
spectrum of energetic particles, that can be under¬ 
stood in terms of an “acceleration zone” near the sheet 
centre, in which the electric field exceeds the magnetic 
field. The escape rate from this zone is then approx¬ 
imated as the time taken by a particle to complete 
one quarter of a revolution around the linking com ¬ 
pone nt of the magnetic field IZenitani and Hoshind 
l299ll| . However, the role of the kink instability ap¬ 
pears to place a relatively modest maximum energy 
limit on the acceleration process. 

5. Summary 

Although the details of the plasma physics remain 
obscure, simple kinematic considerations suggest that 
acceleration at shocks imprints a characteristic power- 
law index on the particle spectrum. In the case of 
relativistic shocks, it is p~ 4 ' 2 , and seems not to be 
sensitive to nonlinear effects, or the effects of magnetic 
field generation at the level of a + ~ %. 

A spectrum consistent with this prediction has been 
identified in a few objects, but observations also show 


that acceleration into a much harder spectrum is 
needed at low energies. Current sheets are in principle 
capable of producing particles with such a spectrum, 
but a full understanding of the way they operate re¬ 
mains a challenging goal. 
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